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By James E. Hubbartt, Henry O. Slone, and Vernon L. Arne 


SUMMARY 


An approximate method for rapid determination of the pressure 
change for subsonic flow of a compressible fluid under the simultaneous 
action of heat transfer, friction, rotation, and area change is de- 
veloped. in the development of this method, the momentum equation was 
approximated and rearranged for a convenient solution employing charts. 
This report presents both the analysis involved in simplifying the mo- 
mentum equation and the charts necessary for obtaining particular solu- 
tions. The charts provide a step-by-step solution which converges to 
an exact solution as the number of steps is increased. An illustrative 
example and comparison with more rigorous numerical solutions with con- 
ditions typical for &ir-cooled turbine blades are included. These com- 
parisons show that the solution converges rapidly to provide кода. 
accuracy. 


INTRODUCTION 


The effective design of ducts to accommodate air flow requires par- 
ticular solutions of the momentum equation for determining the pressure 
changes encountered. Simplified methods which give an accuracy suffi- 
cient for engineering purposes are in demand. An approximate method 
for determining such particular solutions has been developed for the 
one-dimensional flow of a compressible fluid under the influence of heat 
transfer, friction, rotation, and area change. 


A number of studies of one-dimensional flow of a compressible fluid 
have been published. Reference 1 presents a rather complete treatment 
including a form of the momentum equation with heat transfer, friction, 
and area change which is suitable for numerical integration. As an aid 
for expediting the numerical integration, coefficients of the differ- 
entials in the momentum equation of reference 1 have been tabulated as 
"influence" coefficients in reference 2. The analyses of references 1 
and Z have been extended in reference 3 to include rotational forces. 
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Reference 3 also provides a form of the momentun equation suitable for 
numerical integration and tabulates “influence” coefficients for ex- 
pediting the calculations. However, in all cases, the numerical inte- 
grations as employed in references 1 to 5 are tedious and time consum- 
ing. Particular solutions suitable for rapid determination of pressure 
changes for one-dimensional flow with both arbitrary and specific heat- 
input distributions have been presented in references 4 to 7. However, 
these solutions include only heat transfer and friction. In addition, 
reference 8 developed an approximate solution for constant flow area in- 
cluding rotation as well as heat transfer and friction. The approxima- 
tions of reference 8 were based upon experimental data for cooled tur- 
bine blades tested under static conditions. Although the time required 
for a particular solution of the momentum equation and subsequent eval- 
uation of the pressure change was appreciably reduced when the approxi- 
mate solution of reference 8 was utilized, the time reguired for a large 
number of particular solutions, as frequently required for design pro- 
cedures, remains very lengthy. In addition, reference 8 does not account 
for flow-area changes frequently encountered in actual practice. 


Since a method for the rapid determination of the pressure change 
through ducts with heat transfer, friction, rotation, and area change 
is in demand, the momentum equation has been approximated and rearranged 
for a convenient solution employing charts. The purpose of the present 
report is to present the analysis involved in simplifying the momentum 
equation so that a graphical solution is possible and to present the 
charts necessary for obtaining particular solutions for subsonic flow. 
Heat transfer, friction, rotation, and srea changes are all included in 
the solutions. The charts provide a step-by-step solution which con- 
verges to an exact solution as the number of steps is increased and is 
limited only by the graphical accuracy. Comparisons are also made with 
the solutions determined from the more exact numerical integration of 
reference 3 for typical examples of air-cooled turbine blades. 


ANALYSIS 


The general momentum equation for one-dimensional flow is approxi- 
mated and rearranged to form the basis for the construction of charts 
suitable for the determination of the Mach number change with flow under 
the simultaneous action of heat transfer, friction, rotation, and area 
change, A relation between the Mach numbers and pressures is finally 
determined such that charts may be employed to evaluate the pressure 
Changes in a duct once the Mach number changes are known. The energy 
changes may be arbitrarily specified in terms of the total-temperature 


changes and may be evaluated by a procedure such as that given in ref- 
erence 5, 
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Evaluation of Mach Numbers 


The differential form of the general momentum equation for one- 
dimensional flow is (ref. 3) 


aF 
ар - Fgog dx + pV ау + — = 0 (1) 
where x increases in the flow direction (all symbols are defined in 


appendix A). The differential drag force can be expressed by the 
equation 


2 2 2 
_ „ру _ „ру _ Af pV 
dFa = f —— dS = f —— 1 дх = D. 2 A ах (2) 


and the body forces per pound of weight flow, as resulting from rota- 
tion about a central axis, can be expressed by the equation 


2 
or 


If the differential momentum change is replaced by the equivalent ex- 
pression 


oV av = alov?) - v alov) (4) 


and equations (2) to (4) are substituted into equation (1), the momentum 
equation, as integrated between stations x and х+Ах, becomes 


i (жаш (OV) лу 
aP + a(pv) - V d(pv) ~ 
Ру (ov? ).. (pv) 


x 


X+Ax x+tAx 2 
orp dx + a= EY ax = 0 (5) 


Dy 


x x 


The first two terms of equation (5) can be integrated directly. How- 
ever, the last three integrals of equation (5) cannot be evaluated 
directly and therefore are approximated in order to obtain forms which 
can be conveniently solved. The third integral, which exists only in 
the case of area change, and the fifth integral, which expresses the 
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drag forces, may be conveniently approximated by using the integrand as 
the average of the values at the end points. That is, 


(oV). лу 

= (у +V ) 
v alov) ~ XE Гру) адь - (OV) | (6) 
(pV).. 
and 

х+Ах 
ett ы. [time (48) оза G] on 
Ру 2 8 | Dh,x+Ax @ /x+Ax Dn x 2 


x 


In order to facilitate a convenient solution, the fourth integral, 
which expresses the rotational forces in equation (5), will be approxi- 
mated by either  . - 


х+Ах : 
|, arp dx w аёу o Ax (8а) 
х+Ах = | 
1 аётр ах = WT л DA ЛХ (8b) 


depending upon whether the solution is initiated with the conditions at 
the station x ок the station x+\x. Since the integrand of either 
equation (8a) or (8b) is approximated by the.value at a single end point, 
the approximation is in general less accurate than those given by equa- 
tions (6) and (7). Equations (6) to (8) become: identities in the limit 
as Ax approaches zero. 


Because the density and velocity variations from stations x to 
х+Ах are approximated by equations (6), (7), and (8), the energy equa- 
tion must be employed only to determine the total energy change. 


When the approximations given by equations (6), (7), and (8a) are 
substituted into equation (5) and the first two terms of equation (5) 
are integrated, the approximate momentum equation becomes 
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V. + V 
Ру+ду 7 Px + (02) алх - (ОЧ )у - = s [COV алк = (ev) - 
er Ax | xgAx d (ве) tx ew), E 
uc E 2 he B. 2 Jx У TAN 


However, if equation (8b) is used instead of equation (8a), the approxi- 
mate momentum equation becomes 


V + V 
Px+Ax - Ру + (eV2)x+Ax - (oV2)x - аше: [CU - (ву), | - 


4f 2 4f 2 | 
2 Ax х+Ах (m ) x (21 J 
падежа К 2 Aum 2 х+Ах Dn x 2 


If the variables of equations. (9а) and (9b) are changed by employing 
the continuity equation, the relation between the total and static tem- 
peratures, and the equation of state, the two forms of the approximate 
momentum equation can be reduced to the useful relations (appendix B) 


m: Cr Ax 
F (Mews Адк) = Тз (==) coat + — а 949 | 


(1028) 
and 
T or, 
F(M., A.) = 122 EDENA Bins) = — WR] 

(10b) 

where 
Е(М,А) = n 1 - = (a-4) (11) 

M (2 + Y= №) 
(1 ti wey 

G(M) = 2 (12) 


gRM 
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af Ax Ах 
A. = os (= = 1) (13) 


The total-temperature ratio in equations (10) may be evaluated from 
a solution of the energy equation such as presented in reference 3. 
Evaluation of Pressures 


By applying the continuity equation and the equation of state, the 
weight flow can be expressed as 


w = oe VA = ET S MA (15) 
or 
w = — |8 «(^ (16) 
үт" АЕ CUNT 


Since the relation between the total and static temperature ав given 
by the energy equation is 


пе = 1 + 02м. — (17) 


equation (16) can be rewritten in the form 


ї 1/2 
LEE = тё М (a + үс м) / (18) 


Using P'/P = (т'/т)ї/(т-1) to express equation (18) in terms of the 
total pressure gives 


ү+1) 


| aly- 
AT = (te м (1 + Tt yw) к (19) 
P'A R e 
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Equations (18) and (19) are plotted in figure l for ү = 1.40 with 

wv MY PA and w T'/P 'А as the ordinates and M as the abscissa. 
From the Mach number the value of either w T /P ЈА or тутт [PA can 
be determined from figure 1. Then from the given conditions and area, 
P or Р! can be computed. 


CALCULATION PROCEDURE 


Equations (108) and (10b) represent the final forms of the momentum 
equation to be used in solving for the Mach number at either station x 
or x4Ax. If the conditions are known at station x, the Mach number at 
station x+Ax is determined from equation (10a); while if the conditions 
аге known at station x+Ax, the Mach number at station x is determined 
from equation (10b). Increasing values of x, corresponding to а posi- 
tive Ax, must be in the flow direction. A given duct can be divided 
into an arbitrary number of increments in Ax (not necessarily of the 
same length), and the Mach number change through that duct can then be 
computed by progressing stepwise from the end where all conditions are 
known. The size of each increment should depend upon the accuracy de- 
sired. Subsequent examples will indicate the convergence of the solu- 
tion as Ax is decreased. 


Both Е(М, А) and G(M) have been plotted in figure 2 for air with 
Y = 1.40 to expedite the solution. Since ү has a secondary effect 
on the Mach number change and the variation of ү encountered in air 
flow is small, the consideration of other values of y is unwarranted. 


The procedure for determining the Mach number change and thus the 
pressure change with air flowing in a duct is to use equation (10a) or 
(10b) and figures 1 and 2. For example, the procedure for determining 
the Mach number at station х from the conditions at station x4^x is 
to first evaluate Á. ax (eq. (14)). From the value of the Mach number 


at station xx and Á. Ax; the values of F(Mx+Ax,ÁAc-+Ax) and 
G(M., AX) are read from figure 2. Then, with the specified values of 


Or. љуба ду and МТ лт (Ај Али)» FMs A) is calculated 
from equation (10b). From F(M,,A,) and A, (eq. (13)), the Mach num- 


ber at station x is read from figure 2. Once the Mach numbers are 
known, the pressure change can be determined from figure 1. To illus- 
trate this procedure, a typical sample calculation for an air-cooled 
turbine blade including heat transfer, friction, rotation, and area 
change has been carried out in table I with the passage divided into 
two equal parts. Double columns (column 9 and columns 20 to 27) are 
included in the sample calculation setup in table I so that the calcu- 
lations could proceed from either the x or x+Ax station. For each 
of these double colums, the left column is applicable if the solution 
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is initiated at station x+Ax, while the right column is applicable if 

• the solution is initiated at station x. The first 11 columns of table 
I give the conditions which must be prescribed. Although the friction . 
coefficient f is assumed to be the ваше at both the x and х+Ах | 
stations, it may be allowed to vary if greater accuracy seems warranted. 
The remaining columns of table I give the computations required to đe- | 
termine the Mach number at either station x or x4^x. The sample cal- 
culation shown is for the determination of the Mach number at station 

x from conditions at station XMix. 


RESULTS AND DISCUSSION 


In order to evaluate the accuracy of the solution presented herein, 
calculations were made for conditions typical of those occurring for 
air-cooled turbine blades (1) by use of the charts and (2) by the method 
of numerical integration presented in reference 3. The results of these 
calculations employing equation (10b) with 1, 2, and 4 steps with incre- 
ments of equal length are tabulated in table IT and compared with the 
results obtained from the more detailed numerical methods of reference 
3. Three cases are illustrated in table II: (Т) heat transfer and 
friction; (II) heat transfer, friction, and rotation; (III) heat transfer, 
friction, rotation, and area change. For the first two cases, the air 
temperature was assumed to vary linearly from the inlet to the exit. 
However, for the third case, the air temperature was varied exponentially 
to correspond with flow inside a duct heated by crossflow and with con- 
stant inside and outside surface heat-transfer coefficients. The total 
pressures used for evaluating the percentage differences between those 
obtained by the two solutions were computed from figure 1 and the known . 
conditions and areas. 


It is shown in table II that the solutions converge rapidly. In 
addition, the comparison between the two solutions shows that good agree- 
ment is obtained by using only two steps. Even for the third case, the 
computed total pressures agreed within 4 percent for only two steps. In 
shis case the disagreement was essentially eliminated by using four steps. 
Solutions were also made with eight steps, but for the conditions used 
in the examples of table II, the improvement over the use of four steps 
Was probably not within the accuracy of the charts. The approximations 
given by equations (6) to (8) are least accurate at the higher flow Mach 
numbers where slight changes in the total pressure cause large changes 
in the Mach number (fig. 1) and, therefore, large changes in the veloc- 
ity and state conditions. For this reason, the Mach numbers used in the 
examples of table II are intentionally high. It should be pointed out 
that the accuracy of the approximations employed in equations (10) are 
dependent upon Mach number changes, and not the passage length, when ro- 
tational forces are neglected. Thus, the accuracies shown by the results 
of table II are indicative of the accuracies which would be obtained for 
longer passages with similar Mach number changes due to heat transfer, 
friction, and area change alone. | 


CU-2 
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In many cases the air temperature varies exponentially with x as 
in case III of table II. In order to check the effect on the Mach num- 
ber change of using the more easily obtainable linear air temperature 
distribution instead of an exponential variation, the examples presented 
in table II were recomputed as follows: (1) Cases I and II, which were 
calculated for table II with a linear air temperature distribution, were 
determined for an exponential distribution, and (2) case III was com- 
puted for a linear variation instead of the exponential air temperature 
distribution used for table II. The results of these calculations showed 
that the differences between the two variations were negligible and prob- 
ably of the same order of magnitude ав the accuracy of the charts. For 
this reason, it appears that although the air temperature might vary ex- 
ponentially with x, for convenience a linear variation may be assumed 
for evaluating the temperatures needed for the incremental subdivisions. 


CONCLUDING REMARKS 


A method is developed for the rapid determination of the Mach num- 
ber change and pressure change for the flow of a compressible fluid in 
a duct with heat transfer, friction, rotation, and area change. The 
method uses a step-by-step integration which, when applied to specific 
examples typical for air-cooled turbine blades, converged rapidly and 
provided good accuracy with, at most, only four steps. In addition, for 
these specific examples, a negligible difference occurred between using 
& linear or an exponential temperature variation along the duct length. 
This result indicates that the more convenient linear temperature vari- 
ation may be used for evaluating the temperatures for the incremental 
subdivisions if the actual temperature varies exponentially, as fre- 
quently occurs in practice. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 14, 1954 


lO 


NACA 


APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 


_ 1/2 
(1 dira м) i 


flow area, sq ft 
duct length or span, ft 


hydraulic diameter, 44/3, ft 
generalized body force, 10/1Ь of coolant, o^r/g (eq. (3)) 


drag force, 1b 


mE Жш ыры as 


M lot 


friction coefficient 


SEM (eq. (12)) 


standard acceleration due to gravity, 32.2 ft/sec? 
wetted perimeter, ft 

Mach number relative to passage 

static pressure, lb/sq ft abs 

total pressure with respect to passage, 15/sq ft 
gas constant, 55.5 ft-1b/(1b)(?R) 

radius, ft 

surface area, sq ft 

static temperature, OR 


total temperature with respect to passage, ?R 


TN 35150 
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ү velocity relative to passage, ft/sec 
W weight flow, 1b/sec 
x distance measured from passage inlet in flow direction, ft 
Y ratio of specific heats, 1.40 
A iip tas Ax 
x - 2[.—— - 1) (ea. (15)) 
Dx PAX 
AT yx DX Ах 
х +AX 
Naa s - oes - ) (eq. (14) 
D 
h, X-+Ax Ах 
р mass density, slugs/cu ft 
Q a rAx 
mt 
@ angular velocity, radians/sec 
Subscripts: 
е exit of passage 
i inlet of passage 


x distance measured from passage inlet in flow direction, ft 
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DERIVATION OF APPROXIMATE MOMENTUM EQUATION 


Upon collecting terms and employing the relations 


(oV), (>) _ (рү), (OV) scan 


Vx 
(OV) лу (=) = (07) sing tov), МЕЛ “3 a 


and n 


(or), t) -. (ov), (удана Caps) - © ма S 


2 (РҮ) x cx 2 (oV) Ax 


equation (98) becomes 


р P (OV any | 2 2 (27), Fase |, 

aren С i (07) Ax Dp халх 

(ov?) (oV) 4f Ax 
x x+Ax x 2 _ 
4 LI - Or p Ax = 0 (20а) 
2 
and equation (95) becomes 

У Я (оу) л | oo AV AF xt |; 

а & (Vjeme Dh, x+Ax 


(072). V). Af 
int E - 2 es + t - ау Ax Pac ag AX = 0 (205) 


If the continuity equation is employed to obtain the relation 
w = (oVA), = (PVA) aax 


or 


(07) Axa Ax (21) 
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and if the simplifying notations 


£f Ax Ay 
t= EE =) s 
АҒ y AA = J 
= - ——— - 2 14 
ХАХ Dn xx Ко + ( ) 
б: a TAX 
m 
are used, equation (20a) becomes 
ve Ve 
Раље 9 M. (A. Ax ~ 4) pen (A, - 4) - Tipy = 0 
(22a) 
and equation (20b) becomes 
ve ve 
P tig - P. - (гҮ Does - aje VT Ix - &)- АТ л P Ax = О 
(225) 


А more convenient form of equations (22a) and (22b) can be obtained if 
the variables are changed by employing the equation of state, the energy 
equation, and the continuity equation. From the equation of state, 

oV can be expressed as i 


2 


о = a yPM* (23) 


Eum 
gRT 
In addition, the static pressure can be eliminated by using the con- 
tinuity equation and the equation of state as follows: 


м = БРАМ = рт AV (24) 
or 
м ВТ w [ТЕ [Т 1 
POAW =“ № NT M (5, 
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where the total- to static-temperature ratio can be eliminated by using 


the energy equation in the form 


E Ү-1 ма | 
т = 1 + 2 M 
to give the expression 
p.¥ T'R $ 
A 1/2 
d мі + Y- м) 
2 
Finally, application of the continuity equation gives 
гр = eo.» i 
goAV БА V 


mio = 9 D TE 
Ag NrgR N T M 


or, by incorporating equation (17) in equation (28), 


1/2 
` + Y-l м) / 
Ag NTER M 


(17) 


(28) 


(27) 


(28) 


(29) 


If equations (23), (26), and (29) are substituted into equations (22a) 
and (22b) and wi/R/Yg is factored out, the two forms of the approxi- 


mate momentum equation become 


т; ме 
x+Ax 1 | T ду cus - 0) - 


A, 1/2 "7 Uk 40-85 ущ 
"AE 


and 
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15 
Tx 1 | AK (А А, „ә 
== ны eae а ёз 
"e Mese(2 + 152 Мед) 
1/2 
{ту 1 | PE ara) Р TAK (б erga)" Е 
Ах TERT M 4 FAX AL Ay BR Moe Ax 


(30b) 


Upon simplifying, rearranging, and employing the definition of Q, equa- 
tion (30a) becomes 


| от 
rc hate) (а) [round +” о] а 


and equation (30b) becomes 


т! А wor Ах 
ыа [= (252) [rao - HEE ea 
(10b) 
where 

F(M,A) - — —À—;j7l - ne (44) | (11) 

м + r= м2) 

Y-l „2 1/2 

G(M) = u 2 №) (12) 
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TABLE I. - SAMPLE CALCULATION AND SETUP 


pepe pe Теје јаје | e Toa јез | es Tos fon 
Dy x’ 


Ах, | А, | Ахлу, "ax! Th ухх) al Mam, Гул | Уж» | 9, Er "rd = 
Tt | sqft | aq ft ft |+ joe (3) E | e. 
V2) (7)/(8) НЕЕ 


ЕЗГЕ ТҮТТҮ |н jc та |н аза ҮТ ажом 


Мл [Mx | (Мулк) | б(М,„) 
ТАк mi x0 | <o 
(18) (14)272 


zer | awa L9 9 ы ы Бы 


0. | 0.5858 | -0. | -0.5655 / 


DIE 
Q8) ` | 
Ду uawan F(M,; A. F(M, | Ax? уда > 
aer И а [(25) (29)] ane 
х+Ах UT 


JN MEME TE 


(18) 


4f Ax 


OSTS ML VOVN 


LT 
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TABLE IT. - COMPARISON OF RESULTS FROM NUMERICAL AND 
GRAPHICAL SOLUTIONS 
Inlet Mach 


number M4 
from charts 


.1431| .123 


“Heat transfer and friction. 
Heat transfer, friction, and rotation. 


“Heat transfer, friction, rotation, and area change. 
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Figure 1. - Evaluation of welght-flow parameter from Mach number 


(ratio of specific heats, 1.40). (A 152- by 17-in. print of 


this Pig. is attached.) 
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Figure 2. - Evaluation of Е(М,А) and G(M) (ratio of specific 
heats, 1.40). (A 12- by 20-in. print of this fig. is 
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Figure 2. - Continued. Evaluation of F(M,A) and G(M) (ratio of 


specific heats, 1.40). 
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(c) Mach number range, 0.20 to 1.0. 


Figure 2. - Concluded. Evaluation of F(M,A) and G(M) (ratio of specific heats, 1.40). 
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figure 1. - Evaluation of weight-flow parameter from Mach number (ratio of specific beate, 1.40). 
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(a) Mach mumber range, 0.045 to 0.080. 
Figure 2. - Evaluation of F(M,A) and G(M) (ratio of specific heats, 1.40). 
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(b) Mach number range, 0.075 to 0.25. 
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Evaluation of p(M,À) and G(M) (ratio of specific heats, 1.40). 


Figure 2. ~ Continued. 


